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1. Introduction 

The recent interest in phenomenologically viable Technicolor models has driven the explo- 
ration of the phase diagram of SU (N) gauge theories[l]. Investigations have focused on the search 
for asymptotically free theories with infrared fixed points (IRFP) [2], known as the conformal win- 
dow. Such theories have a non-trivial zero in the beta function which leads to conformal behaviour 
in the infrared. 

In particular, Walking Technicolor eases the traditional tension between large standard model 
fermion masses and small flavour changing neutral currents (FCNC) by having the Technicolor 
gauge group run very slowly, walk, between the Technicolor and Extended Technicolor scales. 
This causes an enhancement in {^V^etc, 



where y is the fermion anomalous dimension. Since the standard model fermion masses are pro- 
portional to (*¥ x ¥)etc this allows large fermion masses whilst retaining acceptably small FCNCs. 

Specifically, this requires that the underlying gauge theory be just below the conformal window 
and, equally importantly, that it have a large anomalous dimension (y « 1). Currently, most of the 
effort has gone into investigating theories based on SU (3) with large numbers of fundamental 
fermions or gauge groups with fermions in higher representations [1]. 

Instead, in this work, we study St/ (2) with fundamental fermions. This is attractive as it 
requires fewer fermions to be conformal than SU (3), which leads to a smaller naive S parameter 
and reduced computation time. Theoretical estimates for the lower end of the conformal window 
vary over the range Nf ~ 4 — 8 [3] so we decided to focus on the case of six fermions. We measure 
the running of both the coupling and the fermion mass using the Schrodinger Functional method. 
Unfortunately, if the running is small, it is difficult to distinguish between a truly conformal theory 
and one that is merely running slowly. The running of the mass is easier to determine but requires 
the location of a fixed point to calculate the anomalous dimension. For examples of other theories 
investigated using the Schrodinger Functional see e.g. [4, 5, 6]. 

Subsequent to this talk being given, a fuller account of our findings appeared in [7] and we 
refer the reader to this for more detail. 

2. Formulation 

We define the running coupling g 2 using the Schrodinger Functional (SF) formalism [8, 9]. 
Aside from the fermionic content we follow the method of [5]. The system is defined on a hy- 
percubic lattice of size L with periodic boundary conditions in the spatial directions and Dirichlet 
boundary conditions in the temporal direction. The boundary spatial link matrices are set to: 




(1.1) 



U (x,k)\ t =o = exp[T]Z3a/iL], U (x,k)\ t -L = exp\{% — ri)x^a/iL] 



(2.1) 



with tj = 7r/4[10]. The boundary fermionic fields obey the following relations: 



P + y = 0, yP = aU = 0, 



P-\\r = 0, \jfP+ = at / = L 



(2.2) 
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Figure 1: Data for the running coupling in the Schrodinger Functional Scheme. Simulations are performed 
over a range of bare couplings and lattice sizes L/a. 



with the projection operators defined by P± = ( 1/2) ( 1 ± yo)- 
We define the coupling constant through the relation 



(2.3) 



where k = —24L 2 /a 2 sm(a 2 /L 2 (n — 2r] )) is chosen so that g 2 = g 2 to leading order in perturbation 
theory. This is a non-perturbative definition of the coupling that only depends on L and the lattice 
spacing a. We can then remove the dependence on a by taking the continuum limit. 

We use the Wilson plaquette gauge action, together with fundamental Wilson fermions and an 
RHMC algorithm with four pseudofermions [11]. SF boundary conditions allow us to simulate at 
Kc, which we define as the value of K at which the PCAC mass vanishes. 



mpcAC 



\{dy + d*)f A {L/2) 

2f P (L/2) 



(2.4) 



where /a and fp are the axial and pseudoscalar correlation functions and 5o,^o we t ' ie forwards 
and backwards lattice derivatives respectively. 

As in [12], we use fp to define the pseudoscalar density renormalisation constant: 



Z P (L) = v /3/1//p(L/2) 



(2.5) 



where f\ is a purely boundary term. By observing the scale dependence of Zp we can calculate the 
mass anomalous dimension. However, in order to obtain the correct normalisation this has to be 
measured on lattices with the spatial link matrices at t = and t = L set to unity. 
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Figure 2: Results for o(u)/u using both constant and linear continuum extrapolations. Blue error bars 
indicate statistical error, black error bars include systematic errors. A zero in the beta function is indicated 
by a{u)/u = 1. 



We perform measurements on lattices of size L = 6, 8, 10, 12, 14, 16 for a range of bare cou- 
plings between /3 = 2 and /3 = 8. Measurements to determine K c are made on lattices with unit 
boundary conditions for L < 12 and extrapolated from these for L = 14, 16. Full details of the bare 
parameters used can be found in [7]. 



3. Running Coupling 

As shown in Fig. 1, we measure the SF coupling g 2 for a range of j3,L. It is clear that there is 
little variation with L/a, indicating that the coupling runs very slowly. In order to quantify this we 
calculate the lattice step-scaling function: 

Z(u,s,a/L) =g 2 (go,sL/a)\ g2{g(hL/a)=u . (3.1) 

Then using data from multiple values of L/a we can calculate its continuum limit: 

o(u,s) = lim L(u,s,a/L) (3.2) 

a/L->0 

Throughout this work we set s = 3/2. Using our measured results we interpolate in a/L at each j8 
to find g 2 (go,sL/a) at L = 9l, 10|, 15. The results for L = 6 were found to contain large artifacts 
and are not included in the interpolation or further analysis. For each L we then interpolate in /3 
using the functional form [4] : 

1 _ P 

g 2 (15, L/a) 2N 

See [7] for full details on the optimal parameters and associated errors. 

We are now able to calculate L(u,3/2,a/L) using the fits from Eq. 3.3 for L = 8,9^, 10, 10|. 
We can use these to calculate the continuum value o(u) = ct(m,3/2). We perform both a constant 
continuum extrapolation, using the two values of a/L closest to the continuum limit, and a linear 



4 



SU (2) with Six Fundamental Fermions 



Thomas Pickup 



0.45 



0.40 



N 



0.35 



0.30 



0.25 



0.20 



o 

XX ^ 

n 



■o- 



o 



o L/a=6 




□ L/a=8 




o L/a=10 




a L/a=12 




< L7a=14 




v L7a=16 





"5 



_L 



_L 



_L 



_L 



2.0 3.0 4.0 5.0 

P 



6.0 



7.0 8.0 



Figure 3: Data for the pseudoscalar density renormalisation constant Zp computed in the Schrodinger func- 
tional scheme. 

continuum extrapolation. The results for the constant extrapolation are shown in Fig. 2(a) and for 
the linear extrapolation in Fig. 2(b). The figures show a{u)/u with the blue error bars giving the 
purely statistical error and the black error bars including the systematic errors. 

The eiTors in the linear extrapolation are particularly large and indicate that our data cannot 
provide good constraints on the fits. The results from the constant continuum extrapolation are 
consistent with a fixed point with g 2 > 4.02. However, at our strongest couplings we are still 
unable to rule out the possibility that no fixed point is observed. On the other hand, it is clear that 
at large g 2 a{u) is considerably below the 1-loop prediction. Further data is required in order to 
decrease the errors on the linear extrapolation and thus fully account for the errors in taking the 
continuum limit. 



4. Mass Anomalous Dimension 

We measure the pseudoscalar density renormalisation constant Zp over a a range of j8,L; our 
results are shown in Fig. 3. Unlike the coupling data in Fig. 1, there is noticeable variation with 
L/a at fixed /3, indicating a non-zero anomalous dimension y. 

We extract 7 from the step-scaling function for the mass. We define the lattice version as 



Lp(u,s,a/L) 



Z P (g ,sL/a) 



g 2 (go,L/a)=u 



Zp(go,L/a) 

From this we can define the continuum mass step-scaling function, 

Op(u,s) = lim T,p(u,s,a/L). 

a/L->0 



(4.1) 



(4.2) 
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Figure 4: Results for o>(w)are shown in (a) and the corresponding results for y(u) are shown in (b). The 
blue error bars correspond to the statistical error, the black error bars include all systematic errors and the 
red error bars do not include the continuum extrapolation 

Again, we set s = 3/2 throughout our calculations and discard the L = 6 results. We interpolate in 
L at fixed j3 to obtain results for Zp at L = 9j, 10|, 15. Then for each L we interpolate in j(3 with 
the functional form [5]: 

Z P (P,L/a) = f t c i (j > ) . (4.3) 

Full details on the best fit parameters and associated error can be found in [7]. 

We calculate £p(w, 3/2, a/L) using the fits from Eq. 4.3 in Eq. 4.1 and then perform the con- 
tinuum extrapolation to recover o>(m) = o>(«,3/2). Unlike for the coupling, the errors on both 
linear and constant continuum extrapolation are sufficiently small to allow us to combine them to 
estimate the systematic error due to taking the continuum limit. In Fig. 4(a) we plot G>(w) including 
systematic errors from both the fitting procedure and taking the continuum limit. 

When the coupling is weakly running we can construct an estimator for the anomalous dimen- 
sion, 7, from Op(u) , 

log|<7p(w,j)| 



y{u) 



log |*| 



(4.4) 



This estimator becomes exact at a fixed point [5]. Note also that y itself is only a scheme indepen- 
dent quantity at a fixed point. 

We plot the estimator in Fig. 4(b) using the same error bars as for Op(u). We see that at small 
values of u it closely follows the 1-loop prediction. At stronger coupling both the value and error 
increases and at our strongest couplings becomes compatible with y « 1 as desired for walking 
technicolor. However, at the lower end of our prediction for a fixed point, u = 4.02, 7 could be as 
small as 0.135. 



5. Outlook 

In these proceedings we have presented results for the running of the coupling g 2 and the mass 
anomalous dimension 7 in the Schrodinger Functional scheme for SU (2) with six fundamental 
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We find that our results are consistent with a fixed point at g 2 > 4.02. However, at our strongest 
couplings (g 2 « 5.5) we are unable to rule out that the beta function remains negative and that there 
is no fixed point. In either case, we observe that g 2 runs very slowly across the range investigated. 
Currently, the errors on the linear continuum extrapolation are too large to allow us to calculate our 
systematics due to taking the continuum limit. Further computation is required to address this. 

Determination of the anomalous dimension / is less affected by systematic errors. The primary 
source of error is due to uncertainty over the running of g 2 . Assuming that we do observe a fixed 
point then we find that 0.135 < y < 1.03, with the, phenomenologically interesting, upper range 
only being favoured at the extreme range of the coupling investigated. At weaker couplings our 
results agree well with 1-loop perturbation theory. However, at stronger couplings they favour a 
larger anomalous dimension. 

Currently, the greatest source of uncertainty is due to being unable to perform a linear contin- 
uum extrapolation for g 2 . In order to address this, we are currently performing simulations at larger 
lattice sizes. 
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